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Abstract 



We discuss an outside-inside scenario for a first order quark-hadron transition 
in the Quark-Gluon Plasma (QGP) expected to be produced in heavy ion 
collisions, wherein the entire QGP region itself becomes like a subcritical 
bubble, and starts shrinking. We argue that this shrinking QGP bubble will 
lead to concentration of baryon number in a narrow beam like region in center, 
which can be detected by HBT analysis, or as a raised plateau in the rapidity 
plot of baryon number. 
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Possibility of observing signatures of a quark-hadron transition in relativistic heavy-ion 
collisions has been extensively investigated for sometime now. Lattice calculations allow for 
the possibility that this transition may be of first order. Earlier studies of such a first order 
transition have primarily focused on employing homogeneous nucleation of hadronic bubbles 
in a uniform QGP background. These bubbles expand after nucleation and release latent 
heat which heats up the plasma, suppressing further nucleation [|TJ. 

We discuss an alternative scenario for a first order quark-hadron transition in this paper 
by focusing on the finite size of the QGP region produced in heavy ion collisions. Outside 
this region, one either has a hadronic gas (for example, due to energy density gradient in the 
transverse direction 0), or simply the QCD vacuum in the confining phase. In either case, 
there must be a domain wall at the boundary of the QGP region with a non-zero surface 
tension. For simplicity, we take the region outside this boundary to be (at least a thin layer 
of) the hadronic phase with roughly same temperature as inside. 

We solve the hydrodynamical equations for the longitudinal expansion model of Bjorken 
0, and show that transition can proceed entirely due to collapse of this boundary wall. 
We argue that collapsing boundary wall will lead to concentration of baryons in the central 
region in a narrow beam like region (of thickness 1-2 fm and length of few tens of fm), which 
can be detected by the HBT analysis M of baryons, or by a raised plateau in the rapidity 
plot of baryon number. 

It is important to appreciate that, in contrast to usual QGP bubbles which nucleate by 
thermal fluctuations in a superheated hadronic matter, this entire, bubble like, region of 
QGP does not arise due to any nucleation process. Rather, it forms due to certain initial 
conditions which create a dense, interacting gas of partons in a collision process and this gas 
thermalizes later to form a region of QGP. As thermalization proceeds, the order parameter 
develops a value appropriate for high temperature, leading to the formation of the interface 
at the boundary. 

We start by writing down the hydrodynamical equations governing the evolution of 
energy density e during the longitudinal expansion Q. 

de = (e + p) 
dr r 

Here r denotes the proper time and p is the pressure. [We will use natural system of units 
with h = c = 1.} At temperatures well above the critical temperature, T c , QGP will keep 
expanding in the longitudinal direction without any hindrance. [We neglect any transverse 
expansion.] As the temperature drops due to the longitudinal expansion via eqn.(|I|), the 
outward expansion of the wall (in the local rest frame) will become energetically unfavorable 
at some temperature greater than T c . Let us first estimate this temperature for a simple 
case when the QGP region is spherical with radius R. The free energy associated with such 
a bubble is 

Air 

F = — — [p q (T) — p h (T)]R 3 + ArcR 2 a, (2) 

where p q and ph are the pressures in the QGP and hadronic phases respectively and a is the 
surface tension of the wall. Taking the hadronic gas to consist of three massless pions and 
QGP to consist of gluons and 2 flavors of massless quarks, we have 
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(3) 



where B is the bag constant (which is related to T c by Gibb's condition, p q (T c ) = ph(T c )). 
For T > T c , the critical radius of the QGP bubble is, 

Rc(T) = 2 °" m - (4) 

p q {T)-p h {T) 

A QGP bubble of radius larger than R c will expand while bubbles of smaller sizes will be 
subcritical and collapse. By writing T = T c + AT and assuming AT/T C << 1 (which will 
be the case of interest to us) we can find the expression for AT such that for T < T c + AT 
the radius of the QGP region, R, becomes smaller than the critical radius R c . We find, 

We take a = 50 MeV/fm 2 §. With T c = 170 MeV, we find AT ~ 2MeV, for a spherical 
QGP region with R = 5 fm. So this region starts collapsing when the temperature drops 
below 172 MeV. For T c = 120 MeV, we get AT ~ 5 MeV. Hereafter, we use T c = 170 MeV. 

We now estimate the values of AT suitable for a more realistic cylindrical geometry of 
the QGP region. The free energy of this QGP region (with Rt being the transverse radius 
and Rl the half length) can be written as 

F = -\p q {T) - Ph (T)]27rR 2 T R L + aAnR T R L + 2anR 2 T . (6) 

Here, we have taken the end caps of the cylindrical region to be approximated by flat 
disks of radius Rt- As the dynamics of wall will be governed by local geometry alone, we only 
keep relevant terms in eqn.(^|), for a given portion under consideration, and discard others. 
We first determine the critical size for the middle region, near (longitudinal coordinate ) 
z = 0. For this, we neglect the end portion of the cylinder, i.e. the last term in eqn.(|B|). The 
condition for the transverse expansion to be energetically favored is, 

dF 

<0. (7) 



dR_ 

This gives a critical value of Rt = R c /2 where R c is given by eqn.(|]). Thus for a given 
value of Rt, when the temperature drops below a critical value T c + AT tr then the transverse 
dimension becomes subcritical and the wall should start collapsing. Using arguments as for 
eqn.(|), we find that AT ir = AT sp h/2. 

For the longitudinal motion of the wall, we analyze things in the local frame of the 
plasma. The relevant free energy of a small portion of the cylinder at the end region (with 
length being Az, say) is given by terms in eqn.(^) with R^ replaced by Az. The condition 
for the longitudinal expansion of the wall to be energetically favored is now that the partial 
derivative of F with respect to Az be less than zero (similar to eqn.(|7D). This leads to 
the condition, R T > R c - [Here, it is helpful to think of the end caps to have typical 
curvature determined by Rt, instead of being strictly flat.] Again, for a given value of Rt, 
when the local temperature drops below a critical value T c + ATi ong then the longitudinal 
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dimension becomes subcritical and the wall starts shrinking in that direction. Here we find 
that ATi on g = AT sp h. 

We now come back to eqn.(|l]) for the evolution of the energy density. If the total volume 
of the region (including the QGP region and the hadronic region) is Vo and the volume of 
the QGP region is V q then we can write the average energy density as JT|], 

e(r) = e q ^ + e h (l-^) } (8) 
Vo vo 

where e q and are the energy densities in the QGP phase and the hadronic phase respec- 
tively, 

e = E^T* + B, e h = — T 4 . (9) 
9 30 '10 V ; 

By writing a similar expression, as in eqn.(|8]), for the enthalpy (e + p), eqn.(P leads to 
the evolution equation for the temperature as a function of proper time. Let us first study 
the region near z — 0. For a boost invariant geometry, the same temperature evolution will 
be observed at different values of z. [Again, for central portion, we neglect any possible 
effects of the end cap region.] We consider a cylindrical region of a small length Az around 
z = 0. The outer radius of the cylinder is and the radius of the collapsing wall (bounding 
the QGP region) is denoted by R q . We then have V q = tcR^ A z and V = ■uRq A z. As 
we are neglecting the transverse expansion of the plasma, Rt is time independent while R q 
decreases with the velocity of the wall v (T) moving through plasma. That is, R q = —v(T). 
For v(T) we take the following expression 

v(T) = fl [l - ^] b . (10) 

Here a and b are parameters and T < T c . There are many uncertainties in determining 
the velocity of the wall such as the energy flux through the interface, proper account of 
surface tension etc. ||. We will use two sets of values of these parameters. For one set we 
use a = 3,b = 3/2 (valid for T > 2T c /3, see [[j]||) and for the other set we use a sample 
value, a = 2 and 6=1 (this should be used for T > about 0.7T C , so that v(T) remains less 
than the speed of sound). [Velocities for this second set are larger than those used in the 
literature 0, though linear dependence on T has been discussed in 0. We use this set as 
an example, since we expect that the collapse of the wall will be faster with respect to local 
plasma when one allows for even slightest expansion of the plasma with respect to wall, 
especially near T = T c .} 

There is one subtle point here. The velocity as given by eqn. (|T0|) vanishes for T = T c . 
However, we have argued above that a QGP region of finite size becomes subcritical at a 
temperature larger than T c . Therefore, in this case, the velocity of the wall (with respect to 
the local plasma) should be more accurately represented when we replace T c in eqn.([K]) by 
T c + AT. That is, 
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valid for T < T c + AT, here, AT = AT tr , or ATi ong , depending on whether one is discussing 
the wall motion in the transverse direction or in the longitudinal direction. [We mention 
here that the qualitative aspects of our results do not depend on whether one uses eqn.fllOD 
or eqn.fllTD for wall velocity.] 

Figure 1 shows the evolution of the temperature near z = region. Initial conditions, 
prescribed at (a sample value of) r = 3 fm, are, T(r ) = T c + AT, and R q (ro) = Rt- The 
plots are given only up to those specific values of proper time, when the transverse dimension 
of the QGP region, R q , becomes zero. 

For the evolution of temperature near the end regions of the plasma in the longitudinal 
direction, consider a segment of length Zq(t) at the end portion with the transverse dimension 
of the region being Rt- The location of the end portion of the longitudinally moving wall is 
denoted by z q (r) while the transverse radius of the wall is given by R q {j). The ratio of the 
volume of the QGP region V q to the full (local) region Vq is then given by 

V q _ 7rR q (r) 2 z q (r) 

V ttR%z (t) ■ { ) 

Here, z (t) increases due to plasma expansion with the velocity i = z /t. Wall velocity 
in local rest frame of the plasma is again given by eqn.([T]]), so R q = —v(T) and z q = 
z q /r — v(T) (as both velocities are small). For simplicity, we use a single value of AT 
(= ATi on g) for both directions. With the same initial conditions, as in Fig.l, we then 
determine the evolution of temperature in the end region. We find, that the evolution 
of temperature is almost similar here as for the z = region (Fig.l), apart from minor 
differences, (thus we do not show these plots). For example, corresponding to the dotted 
curve in Fig.l, one finds little larger re- heating for the end portion, while for the thick solid 
curve of Fig.l, one finds smaller reheating. [Note that the temperature obtained from the 
evolution equation represents an average temperature for the cylindrical segment, and hence, 
is valid only for time scales larger than the one given by Rt-] 




FIG. 1. Temperature in the central region as a function of proper time. In both of the figures 
in this paper, thick solid curve and dashed curve correspond to Rt = 7 fm, while thin solid curve 
and dotted curve correspond to Rt = 3 fm. Solid curves correspond to the choice of parameters 
for the wall velocity (eqn.(ll)) as, a = 3,6 = 3/2, while dotted and dashed curves correspond to 
the parameter set a = 2, b = 1. 
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We see from Fig.l that largest reheating temperature is represented by the dotted curve. 
First, T drops sharply to about 0.9T C due to expansion of the plasma. This value is low 
enough to cause nucleation of bubbles in the central region, see 0. Note that the heating, 
and hence suppression of bubble nucleation, will be somewhat larger near the wall as com- 
pared to the center. [We assume that reheating is dominated by the collapsing boundary 
wall due to large surface area, compared to any nucleated bubbles.] For r > about 6 fm, 
the reheating dominates over cooling due to the expansion, and T starts to rise. It reaches 
a value of about 0.98T C by r ~ 16 fm and 0.99T C by r ~ 21 fm. Critical bubble sizes 
corresponding to these temperatures are, R c ~ 3.4 fm and 6.7 fm, respectively (eqn.(|j)). 
The entire size of the region is not very large, or, more importantly, the elapsed time is 
not enough to let any bubbles, nucleated earlier in the interior (when the temperature was 
lower there), to grow to these sizes. [Note that the wall velocity of the bubbles, eqn.(|IO|), 
(or eqn.flTlD with appropriate AT), is very small for temperatures close to T c .} Thus, all 
these bubbles become subcritical and collapse away. The end result being that the phase 
transition to the hadronic phase proceeds entirely due to the collapse of the outer wall. 

For some other parameters, e.g. those corresponding to the plot given by thick solid curve 
in Fig.l, the final reheat temperature is not large enough to completely quench nucleation. 
In this case, the picture of the transition may be somewhat mixed. 

Let us now discuss about some possible observable consequences of this type of scenario 
of phase transition. Here we apply ideas used in the early Universe in the context of baryon 
number inhomogeneity generation at the end of quark-hadron transition ||. Due to large 
masses of baryons in the hadronic phase and much smaller masses of quarks in the QGP 
phase, the baryon number density in the QGP phase is larger than in the hadronic phase. 
The ratio R = n^/n^ of the baryon number densities in the two phases has been estimated 
by many people. A rough estimate may give R ~ 10. Evolution of the value of R in a 
collapsing QGP droplet has also been extensively discussed in literature in the context of 
the early Universe M. It seems reasonable that a good fraction of the net baryon number 
inside the droplet may get trapped inside. We will simply assume this to be the case, at 
this preliminary stage of exploration, instead of getting into details of the mechanism of 
trapping. 

In our model of the phase transition in a heavy ion collision, the shrinking of boundary 
wall leads to a picture which is identical to that of a shrinking QGP droplet in the early 
Universe. Immediate consequence is that the net baryon number of the region (which will 
be expected to be non-zero, though small, even for RHIC and LHC energies) will become 
concentrated near the central region of the plasma due to the collapse of the boundary wall. 
[We mention here that people have discussed formation of strangelets due to strangeness 
enhancement of the QGP droplet in heavy ion collisions by evaporating hadrons |J . How- 
ever, as we show below, the longitudinal collapse of the wall is insignificant by the time 
transition ends, i.e. when transverse collapse of wall is completed. Thus, the picture of 
any, strangeness enriched, exotic lump of matter forming at the center of the collision, or in 
different portions due to local fluctuations, needs to be modified.] 

We calculate the longitudinal collapse of the wall, in the center of mass frame, by tak- 
ing its velocity with respect to local plasma to be given by eqn.(|TTD. In the longitudinal 
expansion model, the plasma velocity at longitudinal distance z is given by z/t 0. Taking 
the end portion of the wall to have initial rapidity equal to some maximum rapidity, the 
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evolution of frame rapidity at wall end is shown in Fig. 2, until the time when transverse 
collapse is completed. The final value of this frame rapidity gives the rapidity interval inside 
which most of the net baryon number may be expected to be concentrated. We see that this 
interval can shrink by up to one unit by the time the transition is completed, i.e. when R q 
becomes zero. 

Rapidity vs proper time 
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FIG. 2. Evolution of frame rapidity at the end portion of the wall due to its longitudinal 
collapse in the center of mass frame. 

Large values of the final rapidity of the wall imply that the longitudinal extent of the 
wall is large when the transverse collapse is completed. This concentrates all (or most) of 
the baryon number in a narrow beam like structure at the center of the cylindrical region. 
Transverse dimension of this region should be of the order of wall thickness, i.e. 1-2 fm, 
while its length can be very large, of the order of r when R q becomes zero in Figs.l and 2. 
Such a structure for baryon source can be easily detected by the HBT analysis. 

Further, the reduction in the frame rapidity at the end portion of wall implies that in the 
rapidity plot of net baryon number, one will expect a raised plateau near zero rapidity with 
width determined by the final rapidity in Fig. 2. Equivalently, one may expect a depression 
in the rapidity plot of antibaryons due to annihilations with increased number of baryons. 
Uniformity of this plateau will depend on whether some (or all) of the fragmentation region 
lies inside the boundary wall initially. 

We conclude by summarizing our main results. We have presented a novel scenario of 
first order quark-hadron phase transition relevant for the situation in heavy ion collisions. 
Here, the phase transition proceeds by the shrinking of the interface surrounding the entire 
QGP region, starting at a temperature slightly above T c itself (due to the QGP region 
becoming subcritical). Collapsing boundary wall heats up the entire region due to the 
release of latent heat which inhibits bubble nucleation in the enclosed region. The reheat 
temperature can be close enough to T c so that if any bubbles were nucleated in the central 
region, they become subcritical and shrink away. Thus, finally, the entire phase transition 
proceeds due to collapse of the boundary wall only. We argue that collapsing boundary wall 
can concentrate the net baryon number of the QGP region in center in a long beam like 
region of very small thickness (of order 1-2 fm), which can be detected by HBT analysis. 
Also, as baryons get concentrated in a smaller rapidity interval, it should be observable as a 
raised plateau in the rapidity plot of baryon number (or a depressed plateau in the rapidity 
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plot of antibaryons). We emphasize that a signature of a pencil shaped source of baryon 
number will not only demonstrate the existence of QGP state, it will also show the transition 
to be of first order, and will provide experimental evidence for the baryon concentration by 
shrinking QGP bubbles as has been proposed for the early Universe. 

Other possible signatures of this picture of phase transition can arise in rapidity plots of 
strange particles, in transverse momentum distributions of particles (e.g. baryons) reflecting 
concentration in the transverse direction, etc. It will be interesting to see if the collapsing 
boundary wall can leave some imprints also on the spectrum of photons and dileptons. 

We would like to thank Pankaj Agrawal, Shashi Phatak, and Supratim Sengupta for 
many useful comments. We are very grateful to Larry McLerran for many useful suggestions, 
especially about the HBT analysis to detect the transverse collapse. 
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